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Abstract. We analyze transport data from a quantum point contact 
(QPC), fabricated on a modulation doped Si/SiGe heterostructurc, to extract 
experimental estimates for the valley splitting. The experimental data are fit to a 
form derived from a valley coupling theory that takes into account the fact that the 
quantum well is grown on a miscut substrate. The results of the fitting analysis are 
compared to results obtained by fitting to a simple phenomenological form; both 
methods indicate that electrostatic and magnetic confinement enhance the valley 
splitting by reducing the lateral spatial extent of the electronic wavefunction. 
Consequently, the valley splitting can be much larger than the spin splitting for 
small magnetic fields. We observe different valley splittings for the two lowest 
orbital modes of the QPC, supporting the notion that when steps are present 
at the quantum well interface, the spatial extent of the wavefunction plays a 
dominant role in determining the valley splitting. 
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1. Introduction 

Quantum computing relies on the existence of a weU defined two-level quantum system, 
or qubit, that is separated in energy from all other nearby levels. There are a number 
of experimental systems where this situation is realized, including NMR [T], optics 
[2] , trapped atoms and ions [5] , and a number of solid state systems [H El E] . Spin 
states in semiconductor quantum dots form a promising qubit candidate because their 
lifetimes can be quite long [71 |H] • Spin measurement and quantum gates have been 
performed in GaAs quantum dot qubits [71IHlini[Tni[ni[Il[Tl[Tl[Til[TSl[I7]. While 
still less developed than their III-V counterparts, significant progress has also occurred 
in silicon quantum dots in recent years [T5 1 fT ^ [ ^ [ ^ [ ^ [ ^ [M l [25 1 [2H l [?7 l l ^ [ ^ 
EOl [311 ISa [Ml El [351136]. 

The advantages of silicon as a host for spin qubits include spin-zero nuclei that are 
naturally abundant and small spin-orbit coupling [BJ [33 [35] , properties that favor long 
spin lifetimes [321 HOI SI] ■ However, silicon also possesses nearly-degenerate conduction 
band minima, or valleys, that could interfere with qubit operations or open additional 
pathways for decoherence [42l [43l [44l [45]. In quantum wells, the valley degeneracy 
is lifted, but the problem of decoherence or interference still exists if the splitting is 
smaller than or comparable to the spin splitting [IS] . 
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Several previous measurements of the valley splitting were made in laterally 
unconfined Si/SiGe two-dimensional electron gases (2DEGs). In these experiments, 
the valley splitting was found to be much smaller than the spin splitting, and to depend 
strongly on the perpendicular magnetic field gll HSl UHl [SOI HIl [SH |33]. Neither 
of these observations can be explained adequately without taking into account the 
presence of atomic steps at the quantum well interface [Ml When steps are 

present, the valley splitting is expected to be larger for wavefunctions that cover fewer 
steps [S31 |S51 [Sni [S71 [S5] . Controlled lateral confinement - for example, in a quantum 
point contact or a quantum dot - therefore provides a means to manipulate the valley 
splitting. 

In this paper, we provide a detailed follow-up to a quantum point contact (QPC) 
experiment reported in j53) . Both this paper and |53) analyze QPC conductivity 
data to obtain the spectrum of energy levels and thus the valley splitting. This 
method for determining valley splitting, which is based on the early work of van 
Wees, requires a formula for the valley splitting as a function of magnetic field |59) . 
In [53) . a simple functional form was used for this dependence. In the present work, 
we analyze the QPC conductivity data using a fitting form for the valley splitting 
that is theoretically derived, using a model for the valley splitting in a QPC geometry 
that includes interfacial steps as well as electrostatic and magnetic confinement. The 
theory yields fitting forms for the valley splitting for each of the magnetic orbitals, 
known as Fock-Darwin (FD) levels. Because multiple FD levels are fit using the same 
parameters, the theoretically justified fitting form requires fewer parameters than the 
phenomenological form used in [53j . 

We compare the results from the different analyses by carrying out a detailed 
statistical analysis to assess the error bars and the range of uncertainty for the fit 
parameters. Even though the asymptotic behaviors of the fitting forms at large 
magnetic fields are different, they yield similar values for the valley splitting in the 
regime where the experimental data are dense. For both methods, the valley splitting 
in a laterally confined quantum device is relatively large (a substantial fraction of 
a millivolt), consistent with the hypothesis that lateral confinement increases the 
valley splitting by reducing the number of interfacial steps covered by the electronic 
wavefunction. 

The paper is organized as follows. In section[2] we review the experimental details 
of the conductance measurements in a QPC. In sectionjSj we discuss the valley splitting 
theory and the fitting forms used in our analysis. In section |4| we present our method 
for extracting the valley splitting. In sections [5] and (Hj we present and discuss the 
results of the fitting analysis. In the appendix, we provide details of the theory of 
valley splitting in a QPC in the presence of a stepped interface. 

2. Experiment and analysis 

2.1. Experiment 

The QPC from [53] was fabricated in a Si/SiGe modulation-doped quantum well grown 
at IBM- Watson in the mid-1990s by ultra-high vacuum chemical vapor deposition on 
a p-type Si substrate at a growth temperature of 560 °C. The method is based on 
the approach reported in [60] . The following growth profile for our sample wafer was 
originally reported in [JS] (referred to therein as wafer MOD42): 1 /im of SiGe, step- 
graded from Si to Sio.75Geo.25; 1 /im of Sio.8Geo.2; 10 nm of Si (the quantum well); 
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Figure 1. (a) Scanning electron micrograph of a top-gated Si/SiGe device, 
operated as a quantum point contact (QPC). In this work the QPC was formed 
between the highlighted gates. The 2-point conductance was measured between 
the source and the drain, (b) Experimentally measured conductance through the 
QPC as a function of gate voltage, for constant magnetic fields of 2 T (solid black), 
4 T (dashed red), and 6 T (dotted blue). 



15 nm of Sio.75Geo.25; 10 nm of Sio. 75660. 25, doped with phosphorus at the density 
Nd = (2-3) X 10i*cm-3; a 4 nm Si cap layer with = (4-10) x lO^Scm-^. The 
resulting two-dimensional electron gas density was 5.7 x 10^^ cm~^ and had a mobility 
of 200,000 cm^/Vs. The nanostructured gates (palladium, deposited by electron beam 
evaporation) were formed on a small mesa etched into the heterostructure. The mesa 
was defined lithographically, and the etching was accomplished with a CF4 reactive 
ion etch. The etch depth was ~ 70 nm. An electron micrograph of the device is shown 
in figure [TJa). The yellow- highlighted palladium Schottky top-gates were used to form 
a quantum point contact by applying negative voltages with respect to the underlying 
electron gas. A source-drain bias of 10 /iV at 112 Hz was applied, and the 2-point 
conductance was measured as the gate voltage was swept from to -1.4 V in steps of 
-0.002 V. A perpendicular magnetic field was stepped from to 8 T in increments of 
0.1 T. 

Figure [l|b) shows typical conductance scans as a function of gate voltage. The 
measured conductance was corrected for a series of background resistances, including 
ohmic contacts, the 2DEG far away from the QPC, and resistors in the filter stage, 
so that the conductance plateaus lined up at multiples of e'^/h. The lowest plateau 
occurs at e^//i, indicating that the valley and spin degeneracies are completely lifted 
in this regime. The observation of a plateau at /h in a Si/SiGe QPC was previously 
reported by Tobben et al [61, and Scappucci et al [SI], and in Si/Si02 inversion layers 
by Wang et al [33] . 

2.2. Spectrocopic method 

Van Wees et al [55] have developed a spectroscopy technique for QPCs, using the fact 
that, at finite magnetic fields, the conductance is quantized into multiples of e^/h. 
Transitions between conductance plateaus occur whenever the minimum energy of 
a one-dimensional subband crosses the Fermi level. Experimentally, the position of 
each transition can be mapped out as a function of gate voltage and magnetic field. 
Theoretically, these transitions can be further mapped onto energy splittings, using a 
model for the QPC subbands. 

An appropriate theoretical expression for the quantum levels in a QPC must 
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take into account spin and orbital effects. In Si, our model must also include valley 
effects, as described in [Appendix A| Treating tfie valley coupling perturbatively, we 
can characterize a single non-degenerate subband by three quantum numbers: the 
transverse confinement mode, or FD level (n), and the spin (s) and valley (v) quantum 
numbers. Here, n can be any non-negative integer, and s and v can take the values 
±f . The energy for subband (n, s, v) is then given by 

En,s,v = ('^ + ^) Vi^^c)^ + {^^0? + (f ) 

1 n'^P 

+ -vA,,{Vg,B) + eVo + —, 

where huJo and hcjc are the electrostatic and magnetic confinement energies, 
respectively, the minimum energy of the confinement potential eVo is determined 
by the gate voltage Vg, gusB is the Zeeman energy, and ti^k'^/2m* is the kinetic 
contribution due to electron motion along the channel. A theoretical expression for 
the valley splitting A„ in orbital n is derived in the appendix. 



2.3. Data fitting window 

Our spectroscopy method relics on the presence of a well-formed QPC, which in turn 
requires full depletion underneath the Schottky gates. For this reason, we use only data 
for which we know a quantum point contact is formed in the channel. For gate voltages 
—0.3 < < V, the conductance is dominated by the quantum Hall effect away from 
the QPC constriction. For gate voltages Vg < —0.3 V there is a clear dependence of the 
conductance on both gate voltage and magnetic field. Therefore we restrict our data 
to the range between -0.3 and -1.4 V, the latter being the most negative gate voltage 
for which data was acquired. We also restrict the data to magnetic fields greater than 
1 T, because the uncertainty in the location of the conductance transitions becomes 



large when i? < 1 T. See sections 4.1 and 4.2 for details 



3. Valleys in silicon/silicon-germanium heterostructures 

Bulk silicon has six degenerate valley minima along the equivalent [001] directions near 
the Brillouin zone boundary at fco ~ 0.82(27r/a), where a is the lattice constant of the 
silicon unit cell. However, silicon grown on a relaxed Si^^Gei-^; virtual substrate is 
strained. This effect partially lifts the sixfold degeneracy by causing the four in-plane 
valleys to rise in energy and the two out-of-plane valleys to fall in energy. The resulting 
splitting is on the order of 200 meV for typical structures [M] . The anisotropic effective 
mass effectively enhances this splitting by 15-20 meV, due to confinement in a quantum 
well, similar to the effect in inversion layers |65| . The lifting of the remaining two- 
fold degeneracy of the z-valleys is the valley splitting of interest in this paper. It is 
caused by an abrupt confinement barrier at the quantum well interface. The resulting 
valley eigenstates involve combinations of fc-states from the two z-valleys. The valley 
splitting arises from the potential energy difference in these two states |56j. 

Early theoretical work predicted a large valley splitting in the range 0.1-1 meV, 
for both Si/SiGe and Si/Si02 heterostructures [66l Ell EH] . These theories assumed a 
perfectly flat quantum well with smooth interfaces, at zero magnetic field, and they 
ignored electron-electron interactions. Later, numerical methods were applied under 
similar assumptions, confirming the previous estimates [IHSSIES]- In the late 1990s, 
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Figure 2. Scaled (dimensionless) valley splitting in a QPC, plotted as a function 
of magnetic field. The plots show the valley splitting for the two lowest FD levels, 
Ao (black) and Ai (blue), (a) Ly = 40 nm, 9 = 0.24°, (b) Ly = 80 nm, 6» = 0.24°, 
(c) Ly = 40 nm, 6 = 0.75°. Note that = 0.24° is consistent with fitting results 
obtained in this paper. For all three panels, we take Lx = 100 nm and ip = Tr/i. 



breakthroughs in the growth of Si/SiGe heterostructures made it possible to study 
new aspects of spin and vaUey physics. Many experiments were performed on Hall 
bars, and reported valley splittings that were much smaller than the predicted values 
[171 13H1 Sni HOI ISIl [SH [S3] . At zero field, most experiments suggested a valley splitting 
at least two orders of magnitude smaller than the theoretical estimates. (Note that an 
anomalously large valley splitting has also been reported for a Si/Si02 quantum well 
[70] . While small enhancements of the valley splitting have previously been attributed 
to many-body effects [7T], the enormous enhancement reported in [70] has not been 
explained theoretically.) 

Such unexpectedly small valley splitting measurements have been attributed to 
the presence of steps at the quantum well interface [S31 [551 ES [57], which cause 
a type of destructive interference that suppresses the valley splitting below its 
theoretical maximum. Electron wave functions that cover many steps experience a 
large suppression in valley splitting. Lateral confinement of the electron wavefunction 
decreases the destructive interference by reducing the number of steps that are covered. 

In a QPC, the electron wavefunction is electrostatically confined by the top-gates 
and is therefore expected to have a larger valley splitting than a Hall bar, even in 
the absence of a magnetic field. Applying a perpendicular magnetic field to the QPC 
should further enhance valley splitting, due to magnetic confinement over the length 
scale 



In Appendix A we develop a theory of valley splitting in a QPC, taking into 



account both electrostatic and magnetic confinement. The theory provides fitting 
forms for the valley splitting in the first two FD levels: 

Ao = f e-^ (3) 
Ai =fe-''|l-25|. (4) 

Here, 



p Aniax 



sm(koLxO cos f) 



k^Lj-O cos ip 



(5) 
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= , (6) 



and 



{ko9 sin ipY 
m*e/h'^ 

where Amax is the theoretical maximum for the valley splitting, obtained at a perfectly 
flat interface, is the length of the QPC, and 6 and ip describe the average tilt 
and orientation of the miscut substrate, respectively. In combination with ([T]), these 
equations provide fitting forms for the QPC conductance transitions. Note that the 
magnetic field dependence in ([s]) and Q predicts an increase in the valley splitting 
with magnetic field, a phenomenon that has been observed repeatedly in Hall bars 
[501 [531 m]. 

This result expresses the valley splitting in a QPC as a function of (i) magnetic 
field, (ii) device parameters, and (iii) heterostructure parameters. Before proceeding 
to the detailed fitting in section |4] it is useful to discuss the main features of ([s]) and 
Q, in order to understand why these forms are physically reasonable and to see how 
some of the basic features can be observed in the data. The full data analysis will 
be discussed in considerable detail in section [4l Here, we simply note that transitions 
between QPC conductance plateaus, such as those in figure [T]jb), occur whenever the 
bottom of an energy subband En^s,v crosses the Fermi level. These transitions can be 
mapped out as a function gate voltage and magnetic field, as shown below in figure [4| 
This transition map contains information about the various terms in ([T]), including 
the valley splitting. 

In Q and Q, the magnetic field and the device parameters appear in the 
valley splitting terms involving ojc and wq, respectively. In ([3|, the valley splitting 
of the lowest FD level increases monotonically with magnetic field. Such behavior 
is manifested in figure |4]ja) because the splitting between conductance transitions 
monotonically increases with field at high fields. The heterostructure parameters 
6 and ip determine the magnitude of the valley splitting in ^ and Q; the valley 
splitting may approach its theoretical maximum A^ax when 9 is small. For typical 
heterostructures, we expect to find A,„ax of order 0.1-1 meV [SS] (recall that the 
Zeeman energy splitting at 2 T is 232 /leV). In figure |4ja), the spacing between 
the conductance transitions at fields in the range 2-4 T is relatively homogeneous, 
indicating that the valley splitting and Zeeman energy splitting in that range are 
comparable. Finally, we note that the exponential functions appearing in ^ and 
Q are physically reasonable. The parameter 6, which appears in the argument, is 
a positive quantity that approaches for large magnetic fields or strong electrostatic 
confinement, causing the valley splitting to saturate. Saturation occurs because, at 
strong enough confinement, further lateral squeezing of the electron wavefunction 
produces diminishing returns in enhancement of the valley splitting. In combination 
with small 9, strong confinement enables the valley splitting to approach its maximum 
value, as we observe. While necessarily general, these observations indicate that ^ 
and Q make some simple predictions that can be understood without detailed fitting. 
To make further progress, a statistical analysis is essential, and we report the results 
of such an approach below. 

In order to make contact with our previous work, we report in parallel an analysis 
of a set of equations for the valley splitting that were used in [53] . These are purely 
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Figure 3. Experimentally measured conductance through a QPC, as a function 
of gate voltage. Results are shown for magnetic fields (a) 3.5 T and (b) 7.4 T. 
Solid, red curves: conductance data. Dashed, black curves: fits to Eq. | [TT| |. 
Dotted, blue curves: fitting residuals. 



phenomenological and are given by 

Ao{V„ B) = ^Al^, + {ABoBr, (8) 

A,{V„B) = ^Al^, + iABiB)^. (9) 

This form was chosen for its simplicity and its linear magnetic field behavior at 
large enough B, with slope Abu, a feature that is consistent with Hall bar data. 
The Avg„ terms, if nonzero, allow for a nonzero valley splitting at i? = 0. In 
section |4| we assume that neither Ay^^ nor Abu depend on magnetic field. All the 
parameters may depend on gate voltage. Because of the phenomenological nature 
of ([s]) and the various terms may not all be statistically relevant. We test 
for this possibility by considering all possible functional variations with the various 
parameters Ay^n and/or Abu set to zero. The statistical relevance of each case can be 
assessed by means of an F-test, which weights the value for the fit by 1/v, where 
p = (number of data points) — (number of fitting parameters) [72] ■ The details of the 
F-test as applied to the parameters in ^ and ^ are presented in [73]. Here, we simply 
note that the most statistically relevant model is the variation where Ay^^ is set to zero 
{i.e., excluded from the model), leaving Ay^o, Abo, and A^i as fitting parameters. In 
the following discussion, all further references to the "phenomenological fitting form" 
refer to this 3-parameter model. 

It is important to note that the theoretical and phenomenological fitting forms 
have different asymptotic behaviors at high magnetic field. At moderate field, in the 
regime where the data are dense (roughly IT < i? < 4 T), we find that the separate 
analyses yield similar values for the valley splitting, as discussed in section [6] Outside 
this field range, extrapolation of either model is uncertain. We emphasize, however, 
that the phenomenological model grows linearly without bound at high fields, and 
is therefore physically unrealistic in such a limit. In contrast, the theoretical model 
provides a physically reasonable saturation behavior. 
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Figure 4. (a) Conductance as a function of gate voltage and magnetic field. 
Solid black lines identify the conductance transitions with [n, s, v) indices (1, 1, 1), 
(1,-1,1), (1,1,-1), (1,-1,-1), (0,1,1), (0,-1,1), and (0,1,-1) (bottom to 
top), (b) Same as (a), but including estimated 1-sigma error bars. For a fixed 
gate voltage {e.g., the dashed line at Vg = —0.5 V), the set of transitions {St a „ } 
defines a system of equations used to extract the valley splitting. 



4. Data fitting analysis 

4.1. Identifying the conductance transitions 

The first step in our spectroscopy analysis is to determine tlie locations of the 
conductance transitions in the QPC as a function of magnetic field and gate voltage. 
A transition between conducting modes occurs whenever an energy minimum of a 
QPC subband crosses the Fermi level Ey- In ([T]), this minimum corresponds to fc = 0. 

We consider a fixed magnetic field. At finite temperature, the conductance 
increases gradually from plateau (A^ — 1) to plateau A, as shown in figurejsj Each step 
can be described by a Fermi-Dirac distribution function [71j , so that the conductance 
is given by 



G = 



N 



1 



h ^ 1 



(10) 



exp[(i;, - EY)/k^T] ■ 

The center of the i"^ step corresponds to a particular gate voltage, such that 
EiiVg) — Ep, where i corresponds to a particular set of indices (n, s, v). The resulting 
transition gate voltage is denoted Vn- 

Since (jlj describes an approximately linear relation between Vg and Ei, the 
conductance data may be fit to 



p N 
i—1 



(11) 



th 



to extract the transition voltages. Here, Wi describes the broadening of the i 
transition in terms of gate voltage. In this way, QPC conductance scans were fit 
for 71 fixed magnetic fields with B > \ T. Typical results are shown in figure [3] For 
further analysis, it is convenient to express the conductance transitions in terms of gate 
voltage, rather than magnetic field. Focusing on the lowest, A < 8, non-degenerate 
subbands, the data were therefore inverted to obtain the transition curves BT.{Vg) 
shown in figure |4] 
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Figure 5. Conductance transitions (solid gray lines) from figure |4] 

Corresponding fits are obtained using (a) a phenomenological fitting form (dashed 
blue line), (b) a theoretical fitting form (dash-dot red line). Note that the 
phenomenological form has a larger number of fitting parameters. (See text.) 



4-2. Transition uncertainty 

If there were no experimental uncertainty, the fitting routine described in section |4.f 
would give the precise location of each conductance transition. In a real system, 
however, there is experimental noise which makes the transition location uncertain. 



Additionally, the simple Fermi-Dirac function used in (111 does not account for 
electron correlations [75], or interference effects as the electron enters or exits the 
QPC [59]. These approximations, which have been left out of our analysis due to 
their complexity, could also, potentially, affect the location of the transition. 

To obtain a quantitative estimate of the uncertainty in Bt^. (Vg), we compute the 
following root-mean-square-error (RMSE) at each conductance step: 



1 

MJ2(Gk-Gntj'. (12) 



Here, M is the number of data points in a single conductance step, and {Gk — Gat^) 
is the fitting residual. For a given transition, the RMSE measures how well the data 
are described by the model we use [72 ], and it provides a quantitative estimate for 
the uncertainty. The outcome encompasses uncertainties due to noise as well as any 
physics not taken into account by the Fermi-Dirac step function. 



Equation ( 12 ) describes the uncertainty in terms of conductance. However, for 
spectroscopy purposes, it is preferable to treat the magnetic field as the independent 
variable, and to express the uncertainty in terms of B. We therefore perform two 
transformations. The first converts conductance uncertainty into voltage uncertainty 
for a given transition location Vt; : 



The second converts voltage uncertainty into magnetic field uncertainty at the 
transition location B^. . Because of the irregular shape of the of the transition curves, 
we use the following discretized prescriptions above and below the transition. 
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Figure 6. Fitting results for the valley splitting parameters obtained from the 
theoretical model, as a function of gate voltage. Shaded regions show the fitting 
uncertainty. 



respectively: 



(14) 
(15) 



Note that the RMSE does not give a precise measurement of the error in our 
model. By comparing the fit to the transition data, however, we find that any 



additional error can be subsumed into (14) and (15) through the factors of 2, which 
have been included. The resulting expressions for provide an upper bound on 

the uncertainty for any physics left out of our model |73) . 



4-3. Energy fitting 

We now perform a second fitting procedure, to extract spectroscopic information from 
the transition maps obtained in the previous section. For a fixed value of Vg in figure[4j 
each conductance transition specifies an equation Ei — Ep, with Ei given in (0. The 
resulting set of (up to) 8 equations is solved simultaneously, obtaining a possible set of 
conductance transitions, Bfit.{Vg), which depend on the fitting parameters eVo, hujQ, 
and any parameters appearing in the valley splitting model. (Recall that A; = 0.) 
In the phenomenological form for the valley splitting, there are 3 fitting parameters, 
AvgOi ^Bo, and A^i. For the theoretical form, there are 2 parameters, T and A. By 
adjusting these parameters, we can minimize the difference between Bfn^ and Bt. by 
minimizing x^, defined as 



N 

E 



-{Bti — -Sfiti)^ 



(16) 



4.4- Error propagation 



The transition uncertainties computed in section 4.2 lead to a range of energy 
uncertainties in section [4.31 To generate error bars for the energy fitting parameters, 
a Monte Carlo error analysis was implemented, in which the fitting procedure was 
repeated many times, using randomly generated instances of the transition fields 
within the uncertainty range. Since the uncertainty is usually dominated by physical 
effects left out of the model, there is no particular reason to believe the uncertainty will 
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Figure 7. Reconstructed valley splitting for the theoretical fitting form, as a 
function of magnetic field and gate voltage, (a) n = FD level, (b) n = 1 FD 
level. 




Magnetic field (T) IVIagnetic field (T) 



Figure 8. Reconstructed valley splitting for the theoretical fitting form. Shaded 
areas indicate the propagated uncertainties, (a) n = FD level, (b) n = I FD 
level. Solid black lines: Vg = -0.5 V. Dashed red lines: Vg = -0.7 V. Dotted 
green lines: Vg = —0.9 V. 



follow a normal distribution. We therefore employed an unweighted distribution to 
generate the transition fields. The uncertainty that we report for our fitting parameters 
corresponds to the standard deviation arising from 50 randomized instances. 



4-5. Best fits 

Figure [5] shows experimental conductance transitions, taken from section 4.1 , together 



with the energy model fits obtained in section 4.3 Results are shown for both the 
phenomenological and theoretical fitting forms. The uncertainties shown here are for 
the transition location, as described in section |4.2| 

Note that energy fitting analysis with the theoretical fitting form was restricted 
to the range —0.914 < < —0.3 V. Outside this range, there are fewer transitions, 
and the equations are under-constrained. The phenomenological fit involves one 
more fitting parameter, so the equations are under-constrained over a wider range. 
Consequently, the procedure is restricted to —0.8 < Vg < —0.3 V. Note that while the 
phenomenological fitting form appears to provide a more consistent match to the data 
in figure [5| it also involves a larger number of fitting parameters. 
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Figure 9. Fitting results for the valley splitting parameters obtained from the 
phenomenological fitting form, as a function of gate voltage. Shaded regions show 
the fitting uncertainty. 



5. Results 

Our results for the theoretical fitting parameters F, A, and Tlluq are shown in figure 6 
Based on these results, we can reconstruct the valley splitting, as shown in figure 7 
The propagated uncertainties for this analysis were discussed in section [4j The 
computed uncertainties are shown in figure [8] for several representative gate voltages. 

As discussed above, no fitting was performed at low fields {B < 1 T). The valley 
splittings plotted in this range are extrapolations, based on the field-independent 
fitting parameters obtained at higher fields. Similarly, the conductance transition 
data are somewhat sparse at high magnetic fields. For the n — level, the valley 
splitting parameters at magnetic fields B > 4.5 T were extrapolated from lower fields. 
For the n = 1 level, parameters were extrapolated in the range B > 2.5 T. 

In figure [9j we show results for the phenomenological fitting parameters, AvgOj 
Abo, and A^i. Figure [T0| shows the reconstructed valley splitting. In figure pT] we 
show the propagated uncertainty for the valley splitting at several fixed gate voltages. 
In figures [TO] and [TT] we also include low and high-field extrapolations of our fitting 
results, as we did for the theory plots. 



The valley splitting results shown in figures 10 and 11 agree fairly well with those 
obtained in Goswami et al. [53]. We attribute the slight differences to the Fermi- 
Dirac fitting procedure, used here to determine the conductance transitions, and to 
our restriction of the energy fitting to the range i? > 1 T. In the voltage range 
—0.7 < < —0.3 V, the same number of transitions were fit in both analyses and 
the reconstructed valley splittings were found to differ by less than 0.3 meV. For 
Vg < —0.7 V, the valley splittings reported here differ by up to about 1.0 meV from 
those in Goswami et al, because a different number of transitions were fit. 



6. Discussion 



6.1. Comparison of fitting results 

In the context of quantum computing, we are particularly interested in the overall 
magnitude of the valley splitting. It is therefore essential to know whether the results 
of our analysis depend sensitively on the fitting form that we use. In figure |12[ we 
compare the reconstructed valley splittings for the phenomenological and theoretical 
fitting forms. In the region where the transition data are dense, corresponding to 
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Figure 10. Reconstructed valley splitting for the phenomenological fitting form, 
as a function of magnetic field and gate voltage, (a) n = FD level, (b) n = 1 
FD level. 




1 2 3 4 5 6 7 
Magnetic field (T) 



1 2 3 4 5 6 7 
IVIagnetic field (T) 



Figure 11. Reconstructed valley splitting for the phenomenological fitting form. 
Shaded areas indicate the propagated uncertainties, (a) n = FD level, (b) n = 1 
FD level. Solid black fines: Vg = -0.5 V. Dashed red lines: Vg = -0.7 V. Dotted 
green lines: = —0.9 V. 



1 < i? < 4.5 T for n = and 1 < i? < 2 T for rt = 1, we observe that the results 
do not depend sensitively on the fitting forms, and the two models yield results that 
differ by no more than 0.5 meV. This is an important conclusion for our analysis: the 
extracted valley splitting is not strongly model dependent in the regime where the 
data are dense. 



For high magnetic fields, the reconstructed valley splittings in figure 12 
exhibit very different behaviors. The theoretical fitting form saturates while 
the phenomenological fitting form increases linearly. We have argued that the 
extrapolation of the phenomenological form must be incorrect at high field, on physical 
grounds. The theoretical form does not suffer from this problem. 

It is interesting to note that the low-field extrapolations obtained from the two 
fitting forms agree in their main features: a significant zero- field valley splitting for the 
n = level, and a relatively small splitting for the n = 1 level. The non-monotonic 
behavior of the n = 1 FD level is a conspicuous theoretical feature, and the same 
feature was observed in a previous numerical analysis '76'. In the present experiment, 
such non-monotonic behavior occurs in the field range i3 < 1 T, where the transitions 
are not well-resolved, precluding the possibility of confirming such behavior here. 
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Figure 12. Valley splitting results for the phenomenological fitting form (dashed 
blue) compared to the theoretical fitting form (dash-dot red), for three gate 
voltages. The top row shows the valley splitting for the n = FD level at (a) 
Vg = -0.5 V, (b) Vg = -0.7 V, and (c) Vg = -0.9 V. The bottom row shows the 
valley splitting for the n = 1 FD level at (d) Vg = -0.5 V, (e) Vg = -0.7 V, and 
(f) Vg = —0.9 V. Shaded areas indicate propagated uncertainty. 



6.2. Characterizing the quantum well interface 

The theoretical expressions for the valley splitting, obtained in the appendix, 
incorporate materials parameters describing the substrate miscut angle 9 with respect 
to [001], and the relative orientation of the QPC device (p with respect to [100]. The 
relative orientation of the QPC to the atomic steps at the interface could not be 
determined for this sample; for similar samples, the average orientation was about 
45°. By using our transition fitting procedure with the theoretical fitting forms ([s]) 
and (|4]), it is possible to place some bounds on the values of 6 and ip for this sample 
from the experimental data. 

The main obstacle for determining 6 and (p lies in ([5|, which contains an additional 
device parameter, Amax, the valley splitting 'maximum' for a perfectly flat quantum 
well interface. Amax can be estimated theoretically [56] using the areal density of our 
2DEG (n ~ 5.72 x 10" cm"^), giving A^ax - 0.5 meV. We can then estimate f 
as follows. In we note that sin(a;)/a; has a maximum value of 1. f should then 
have a maximum value of 0.3 meV. Our experimental estimates for F, however, range 
between 0.9 and 1.8 meV, as shown in figure |6] This discrepancy between experiment 
and theory arises in other 2DEG valley splitting experiments [70l [71] . The effect is 
usually attributed to many-body correlations, which are not present in single-electron 
theories. Indeed, our results are consistent with other recent observations of enhanced 
valley sphtting |i50j . 

Since the value of Amax is not well known, it becomes impossible to determine 
6 and ip independently. However, we can place an approximate bound on 6, using 



Valley splitting in a Si/SiGe QPC 



15 




cp (degrees) Gate voltage (V) 



Figure 13. (a) Tilt angle 6 as a function of orientation angle ifi, for the 
case Vg = —0.5 V. (b) Tilt angle as a function of gate voltage, for the case 
sin^((/p) = 1/2. Shaded areas indicate uncertainty. 



i7\ and the fitting results for A. For example, the dependence of 6 on f is shown in 



figure 13 With the exception of the very smallest rotation angles {tp < 20°), we see 



that 9 varies slowly in the range 0.15 < 6* < 0.45°. Thus, a characteristic value of 9 



is obtained by setting sin ip — 1/2 in ([7|). The result, shown in figure 13 ^b), suggests 
that ~ 0.25°. 

6.3. Dependence of valley splitting on the Fock-Darwin level 

A main conclusion of this paper, and of several recent papers (531 ISSl liSl [S71 [55] , 
is that in the presence of a miscut substrate or a rough quantum well interface, the 
observed valley splitting will depend sensitively on the lateral extent of the electron 
wavefunction. In a QPC, the lateral extent of the wavefunction depends on both the 
electrostatic confinement due to the top-gates, and the magnetic confinement from 
a perpendicular i?-field. In figure [12] such behavior can be observed clearly for the 
n = FD level as a function of magnetic field. For higher magnetic fields, the magnetic 
confinement increases and valley splitting is enhanced. This trend does not depend on 
which fitting form we use in our analysis. At low fields, the valley splitting is non-zero, 
due to electrostatic confinement. 



On the other hand, by comparing panels (a)-(c) in figure 12 we note that both 
fitting forms agree that the valley splitting decreases with increasing electrostatic 
confinement. This is an experimental fact, which cannot been explained by the theory 
presented in the appendix. Here, we attribute this behavior to electron correlation 
effects. For weak electrostatic confinement, the electron density in the QPC is 
relatively high, since there are many electron modes with A; > lying below the 
Fermi level. As argued elsewhere, such electron correlations should enhance the valley 
splitting [501 [73 ITU- 

The behavior of the n — 1 FD level is generally consistent with these confinement 
arguments. In figure [12[ similar trends are observed for both n = and n — 1. 
However, at low fields, the valley splitting behavior is complicated by the nontrivial 
nodal structure of the n = 1 wavef unctions, as apparent in the theoretical fitting 
form. Generally, one might expect n — 1 wavefunctions to experience a smaller valley 
splitting than n = 0, since the wavefunctions are more extended. Such behavior is 
observed in the theoretical fitting form at high fields. At low fields, however, the nodal 
structure of the wavefunctions interferes with our expectations. In fact, the low-field 
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valley splitting in panel (f) is actually larger than panel (c). It is worth noting that the 
valley splitting for the n = 1 FD level can actually be tuned to zero by the application 
of a magnetic field. 



7. Conclusions 



In this paper, we extracted the valley splitting from conductance measurements in a 
Si/SiGe QPC using two different fitting forms: a phenomenological form similar to the 
one used in 'S3' and a theoretical form developed here. For both models, the results 
support the hypothesis that steps at the quantum well interface lead to a dependence 
of the valley splitting on the spatial extent of the electron wavefunction. For the 
lowest transverse mode of the QPC, the valley splitting ranges from about 0.5 meV 
at zero magnetic field to 1.5 meV at high fields. For the first excited mode, the valley 
splitting varies non-monotonically and can be tuned from to about 1 meV. 
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Appendix A. Theory of valley splitting in a quantum point contact 



In this appendix, we develop the valley splitting theory for a single electron confined in 
QPC geometry, in the presence of a tilted quantum well. The effective mass envelope 
function equation for the confined electron is 

" 3 / ^ s 2 



E- 



2m„ 



'ih- 



d_ 

'dr. 



(A.i) 



+Vc^w{z') + V{x',y') 



F{r') = EF{r'). 



Here, 



in figure Al 



the primed coordinate system refers to the plane of the interface, as shown 
We have made the usual assumption that the vertical and lateral 



confinement potentials \Vqsn{z') and V{x',y'), respectively] are separable. 

The transverse and lateral effective masses are given by mi ~ = ™t = 0-19 toq 
and TO3 = m\ = 0.92 toq, respectively. Note that the tilted interface induces off- 
diagonal terms in the inverse effective mass tensor, in the primed coordinate system. 
However, for typical tilt angles, these terms comprise less than 3% of the diagonal 
elements, and can be treated perturbatively |56j . In the present work, we ignore 
them. 

Since the atomic steps are much smaller than any effective mass length scale, we 
can make the approximation that the quantum well potential Vqw(-2') is smoothly 
tilted (i.e., not step-like). In the perpendicular field geometry, B = Bz' , separation 




Figure Al. Schematic QPC geometry, including steps, (a) Top view of QPC. 
Lx and Ly are the length and width of the QPC, respectively. Interfacial steps are 
shown as diagonal lines, rotated from the channel axis by angle The downhill 
direction, perpendicular to the uniform steps, corresponds to p' . (b) Cross-section 
view of the steps, with miscut angle 9. Heterostructure layers, from bottom to top, 
correspond to substrate, quantum well, and spacer. Doping layer and top-gates 
are not pictured. 



of variables leads to solutions of form F{r') = Fxy{x' ,y')Fz{z'), where Fz{z') is the 
quantum well wavefunction and F,j.y{x',y') is the lateral wavefunction. 

We consider the specific QPC geometry shown in figure |A1| The confinement 
along y' is modelled as a parabolic potential. In the x' direction, there is no 
confinement, but we assume periodic boundary conditions over the distance L^- The 
step density is taken to be uniform, with step orientation given by the angle (f relative 
to the QPC geometry. 

We now solve (A.ll by what amounts to a Fock-Darwin (FD) analysis [771 178) . 
The zero-field channel width Ly is related to the oscillator energy hujQ as 



There is an analogous relation between the cyclotron frequency, ujc — e\B\/mt, and 
the magnetic length: 

- - si- 

Working in the primed coordinate system, we then have 

H=-^ {-ihV + eAf + \m,u:ly'^ + Vq^{z'). (A.4) 

The quantum well is perpendicular to and we define the position of the top interface 
to be at z' = 0. As typical in a 2DEG, we assume a strong modulation doping field so 
that the electron wavefunction feels only the top interface of the quantum well. The 
resulting valley coupling potential is given by [S5] 

V^{-J)^v^6{z'), (A.5) 

which we shall treat via perturbation theory. The coupling strength depends on the 
conduction band offset at the quantum well, Ai?c, as ~ 5.4 x 10^^^Ai?c. Here 
is in units of eV-m when Ai?c is in units of eV. Note that we have used the more 
accurate estimate of from |55] , based on the many-band tight binding model, rather 
than the two-band model. 
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As noted above, the problem is separable in the (x',y') and z' variables. We now 
solve the {x',y') problem. For the present geometry, with the magnetic field oriented 
perpendicularly to the interface, it is convenient to use the Landau gauge for the vector 
potential: A = —By' x' . Except where noted, we henceforth drop the prime notation. 
The resulting Hamiltonian, involving the lateral variables, is given by 



H = -- 



2mt 



d 



1 



(A.6) 



In (A.6 1, our choice of gauge permits a further separation of variables. We consider 



the ansatz F{x,y) — e:>q3{ikx)u{y) , where we have dropped the xy subscript. The 
eigenvalue equation for u{y) becomes 



^ 

2mt + ojg 2TOt dy^ 



(A.7) 



where 



Vk 



In the X direction, the periodic boundary condition leads to the quantization 

Lx 

where j is an integer. 



Solutions of (A.7) are given by 



Ekn 



2to(- lu 



%+(^+]^nuj (n = 0,l,...), 



and 



where 



Fknix,y) 



2 2 

UJ = UJr. 



exp{ikx) 
y/nLxL2''n 



(A.8) 
(A.9) 

(A.IO) 
(A.ll) 



'■B ' 

and Hn {x) is a Hermite polynomial TS] . The solutions described above are known as 
Fock-Darwin levels. 

The valley splitting is now computed from the prescription [56] |f] 



dr'' exp(-2ifcoz) |Ffe„(r', cj)')F,iz')f K(r') 



(A.12) 



where we have reverted to primed coordinates since the exponential factor exp(— 2ifcoz) 
(i.e., the valley phase) is defined with respect to the crystallographic (unprimed) axes. 
To complete the integral, we need to express the valley phase in terms of primed 
coordinates. In the continuum approximation (i.e., not treating the steps as discrete), 
the following relation applies to the interface: 

z ~ z' + 9y' simp — 9x' cos (p. (A. 13) 



f Absolute values have been used inside the integrand, to account for complex wavefunctions. 



Valley splitting in a Si/SiGe QPC 



19 



Applying (A.5), we obtain 



dx / dy (A. 14) 



X exp[2i/co^'(a;cos(/7 - ysinip)] \Fknix,y)\^ 



where we have again dropped the prime notation. 

In (A. 14 1, the parameter A^ax — 2uv^z(0) corresponds to the theoretical 
maximum for the valley splitting in a perfectly flat quantum well. An estimate for 
Aniax requires knowledge of the quantum well wavefunction Fz{z). The wavefunction 
can be computed analytically for a 2DEG in the Hartree approximation, as described 
in |56j . As usual, we adopt the lowest subband approximation. Expressing the 2DEG 
density n in units of lO^-'^cm"^, we obtain Amax — 0.086n, in units of meV. However, as 



described in section 6.2 this procedure leaves out certain many-body effects, leading 
to an underestimate of the valley splitting. In the following analysis, we circumvent 
this problem by folding Amax into the valley splitting normalization, since it appears 
as a prefactor in all our valley splitting expressions. 

Defining the scaled valley splitting as En = An/Amax, we can obtain results for 
several FD levels: 



El = r(6i,^)e-''|l - 2b\, 

£■2 = T{e, ip)e-^\l - (2 + V2)b\ |1 - (2 - V2)b\, 
where we have defined 



1 



sm{koLx0 cos ip) 



koL^O cos ip 



and 



{koLO sin(^)^ 



(A.15) 
(A.16) 
(A.17) 

(A.18) 
(A.19) 



The geometrical factor r{6, (p) oscillates rapidly as a function of 6 and (p, reflecting 
the fact that 9 and (p control the number of steps inside the channel of the QPC. For 
the fitting analysis presented in the main text, the equations derived above can be 
reduced to simple fitting forms, as given in (|3f and @ 



A typical picture of valley splitting in a QPC is obtained by evaluating (A.15)- 
(A.17) for several realistic geometries, as shown above in figure[2] Note that only the 
n = 0, 1 FD levels are plotted. Comparing the three panels, we observe that the tilt 
angle i? plays a strong role in determining the magnetic field scale, since it is directly 
related to the number of steps covered by the wavefunction. The width of the QPC 
plays a strong role at low fields, where the confinement is predominantly electrostatic, 
rather than magnetic. At higher magnetic fields, we observe oscillations in the valley 
splitting caused by the nodal structure of the wavefunctions when ti > 0. These 
oscillations were first reported in [76] . Indeed, when we calculate the valley splitting 
for the device parameters used in |76| , we obtain a very similar magnitude and field 
dependence. However, the latter work does not take into account the electrostatic 
confinement of the QPC. 
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